Dark-Matter Admixed Neutron Stars 
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We study the hydrostatic equilibrium configuration of an admixture of degenerate dark matter 
and normal nuclear matter by using a general relativistic two-fluid formalism. We consider non- 
self-annihilating dark matter particles of mass ~ 1 GeV. The mass-radius relations and moments of 
inertia of these dark-matter admixed neutron stars are investigated and the stability of these stars 
is demonstrated by performing a radial perturbation analysis. We flnd a new class of compact stars 
which consists of a small normal matter core with radius of a few km embedded in a ten-kilometer- 
sized dark matter halo. These stellar objects may be observed as extraordinarily small neutron stars 
that are incompatible with realistic nuclear matter models. 

PACS numbers: 95.35.-|-d, 97.60.Jd, 



Introduction. — By now the existence of dark matter 
(DM) has been well established, with a large amount of 
evidences such as galactic rotation curves, cosmological 
structure and gravitational lensing. However, the proper- 
ties of DM including their mass and interactions are still 
unknown. It is thus of great interest to constrain the 
properties of DM through direct or indirect methods. 

Direct methods search for the signals of DM-nucleus 
scatterings in Earth-based detectors. The latest exper- 
imental results are not conclusive. The data from the 
DAMA [H and CoGeNT experiments are consistent 
with detecting light DM particles with mass ^ 10 GeV, 
which are incompatible with the null results from CDMS 
[i| and XENON Nevertheless, it has recently been 
suggested that isospin-violating DM may be the key to 
reconciling the experimental results Q . On the other 
hand, indirect methods are based on the effects of DM 
on the properties of stellar objects such as the Sun. For 
example, the effects of low-mass (^ 5 GeV) asymmetric 
DM particles on the solar composition, oscillations, and 
neutrino fluxes have been considered recently 

One indirect method that is gaining attention in recent 
years is to study the effects of DM on compact stars. The 
effects due to different DM models have been considered. 
For example, self-annihilating DM inside compact stars 
can heat the stars, and hence affect the cooling properties 
of compact stars [l3, [HI ■ On the other hand, non-self- 
annihilating DM, such as asymmetric DM [l^ and mirror 
matter , would simply accumulate inside the stars and 
affect the stellar structure. Constraints have been set by 
connecting the observed properties of compact stars with 
DM parameters. 

It should be noted that neutron stars with a DM 
core are inherently two-fluid systems where the normal 
matter (NM) and DM couple essentially only through 
gravity. The technique used in recent studies of the 
structure of these dark-matter admixed neutron stars 
(DANS) is based on an ad hoc separation of the Tolman- 
Oppenheimer-Volkoff (TOV) equation into two different 
sets for the normal and dark components inside the star 
[Tsl [T3 [. This approach is motivated by the similarity 



of the structure equations between the relativistic and 
Newtonian ones, but it is not derived from first principle. 
In fact, a general relativistic two-fluid formalism is avail- 
able [l5l| and has been employed in the study of superfluid 
neutron stars (e.g., ,16, 17]), where the two fluids are nor- 
mal and superfluid nuclear matter. This approach is not 
only more appealing from a theoretical point of view but 
also able to extend easily the study of dynamical proper- 
ties of these stars in a self-consistent general relativistic 
framework. Here we study the structure and stability of 
DANS in general. Besides the scenario where a DM core 
exists inside a neutron star, we also study the scenario 
where NM is in the core of a DM dominated compact 
star. It should be pointed out that the main focus of this 
paper is to study the equilibrium properties and obser- 
vational signatures of these theoretical objects (see also 
[l8| for a study of compact stars made of fermionic DM). 
The formation process of these objects requires further 
investigation. 



Formulation. — To study a two-fluid compact star, we 
adopt the formulation given in [l6| . which was initially 
constructed to study general relativistic superfluid neu- 
tron stars. Here we shall briefly summarize the formal- 
ism and refer the reader to 16] for more details. The 
central quantity of the two-fluid formalism is the master 
function A(n^,p^, x^), which is formed by three scalars, 
— —rian"', = —pap", and ~ —riap"'. The four 
vectors n" and p" are the conserved NM and DM num- 
ber density currents respectively. The master function is 
a two-fluid analog of the equation of state (EOS) and —A 
is taken to be the thermodynamic energy density. 



For a static and spherically symmetric spacetime ds'^ = 
-e^'^^'^dt^ + e^^^'Ur'^ + r'^{de^ + sin^ Odcj)^), the structure 
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FIG. 1: Mass-radius relations for different amount of DM 
specified by the parameter e (see text). The DM particle 
mass is mx = 1 GeV. The density profiles for the ordinary 
neutron star (e = 0) and DM dominated star (e — 0.8) at the 
point X are shown in Fig. [2l 
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FIG. 2: Upper panel: Density profile for the ordinary neutron 
star (e = 0) at the point X in Fig. [T] Lower panel: NM and 
DM density profiles for the DM dominated star (e — 0.8) at 
the same point. 



equations for a two-fluid compact star are given by |16| 
Ay + By + ^{Bn + Apy = 0, 

CoV + A"on' + ^{An + Cp^ = 0, 

1 - ^ 

A = Svrre^A, 

r 

, 1 — 

f = — + 87rre^^', (1) 

r 

where the primes indicate derivative with respect to r, 
and the coeSicients A, B, C, Aq, Bq, and Cq are func- 
tions of the master function. Their expressions are given 
by Eqs. (3) and (25) in [l^. The generalized pressure ^ 
is computed by Eq. (18) in ^1^. 

The EOS information P = P{p) (with P and p be- 
ing the pressure and energy density respectively) needed 
in the standard relativistic-star calculation based on the 
TOV equation is now replaced by the master function 
A(n^,p^, a;^). We assume no interaction between NM and 
DM except for gravitation. The master function does not 
depend on the scalar = —Uap" and is separable in the 
sense that A(n^,p2) = Anm("^) + Adm(p^), ^Nuin'^) 
and Adm(p^) being the negative of energy densities of 
NM and DM at a given number density respectively. We 
use the APR EOS 1^ for NM and assume that the DM 
component of the star is formed by non-self-annihilating 
DM governed by an ideal Fermi gas. As discussed ear- 
lier, DM candidates in the mass range of a few GeV are of 
great interest recently. We shall thus consider fermionic 
DM particles of mass mx ^ 1 GeV in this work. 

Results. — In Fig. [T]we show the mass-radius relations 
of DANS for different amount of DM specified by the 
parameter e — Mdm/(-Mnm -I- My)m)- Here Mnm and 
Mdm are calculated by the product of the particle mass 
and total number of particles for NM and DM in the 
star. They may be referred to as the baryonic masses 
for NM and DM (though it should be noted that DM 



is nonbaryonic). In the figure, M is the gravitational 
mass and R is the radius of the star. The DM particle 
mass mx = 1 GeV is fixed. The case e = corresponds 
to ordinary neutron star models constructed using the 
APR EOS without DM. We see that the existence of a 
DM core would lower the maximum stable mass allowed 
by this EOS. The DANS would also have smaller radii. 
For the maximum stable mass configuration, M and R 
are decreased by about 35% and 9% respectively as e 
changes from to 0.2. For the case e = 0.8, the stars are 
DM dominated compact stars, with NM concentrated in 
the core. It is seen clearly that the mass-radius relation 
of these stars is different from that of NM dominated 
stars qualitatively. 

In Fig. [1] the circle marked by X is the intersection 
point between the curves e = and 0.8. While the or- 
dinary neutron star (e = 0) and dark-matter dominated 
star (e = 0.8) at X have the same M and R, their in- 
ternal structures in fact differ significantly. Fig. [5] shows 
the density profiles of these two stellar configurations, 
the upper (lower) panel corresponding to the ordinary 
neutron star (DANS) model. For the case e = 0.8 (lower 
panel), it is seen that a small NM core is embedded in a 
ten-kilometer sized DM halo. Since M and R of the two 
stars are the same, it would seem impossible to distin- 
guish them based on their gravitational effects on other 
nearby stellar objects. However, the visible radius of the 
DM dominated star (defined by the radius of the NM 
core) is 0.56i? and the total mass enclosed in the NM core 
is 0.72M. The two stars can be distinguished by measur- 
ing the gravitational redshift of spectral lines, since that 
produced near the surface of the NM core will be about 
30% larger than that of ordinary neutron stars. 

In order to check the stability of these stars, we have 
solved the set of equations for radial perturbations of a 
two-fluid compact stars developed in [1^ . Similar to the 
one-fluid case, the problem is to solve for the eigenvalues 
o;^, where w is the oscillation frequency of the star. We 
will present the details of our calculations and analysis 
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FIG. 3: Gravitational mass M (upper panel) and squared 
frequency of the fundamental mode ojq (lower panel) are plot- 
ted against the central energy density for the DM dominated 
sequence (e — 0.8) in Fig. [T] 



elsewhere. Here we show our main results in Fig.[3l where 
M is plotted against the central energy density in the up- 
per panel for the DM dominated sequence (e = 0.8) in 
Fig.[T] The lower panel plots the squared frequency of the 
fundamental mode Wq for the same sequence. Similar to 
the one-fluid study for the stability of ordinary neutron 
stars, the point Wg = marks the onset of instability. 
Fig. [3] shows that ojq passes through zero at the central 
density corresponding to the maximum mass configura- 
tion. Beyond this critical central density, the stars are 
unstable against radial perturbations. For lower central 
densities, such as the DM dominated star at the point X 
in Fig. [1] the stellar configurations are all stable. 

Besides the gravitational mass and radius, it is also in- 
teresting to consider the moment of inertia / of DANS 
since it is measurable and plays an important role in the 
physics of neutron stars. Bejger and Haensel [13] discov- 
ered an (approximately) EOS- independent formula relat- 
ing /, M and R. For ordinary neutron stars, they found 
that 



z/(0.1 + 2z) 
2(l + 5z)/9 



if z < 0.1, 
if z > 0.1, 



(2) 



where the scaled moment of inertia / = I /MR^ and 
z = (M /MQ){km/ R). This universal formula was ob- 
tained by fitting a large set of realistic EOS models for nu- 
clear matter, including the APR EOS used in this work. 
The moment of inertia of a rotating star in general rela- 
tivity is commonly defined by / = J/il, where J and Q, 
are the angular momentum and angular velocity respec- 
tively. In the slow rotation limit, J scales with f2 linearly 
and / is determined only by the nonrotating background 
quantities [2l| . 

In the two-fluid case, individual NM and DM angu- 
lar momenta (Jnm and Jdm) can be defined [17[. In 
the slow rotation limit, the moments of inertia can also 
be defined by /nm = -'^nm/^^nm and /dm = Jdm/^^dm, 
where JInm and r^oM are, respectively, the angular ve- 
locity of NM and DM. /nm and /dm depend on the 
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FIG. 4: Scaled moment of inertia I/MB? is plotted against 
compactness M/R for the DANS sequences shown in Fig. [T] 
Each sequence is labeled by the parameter e. The dashed line 
represents Eq. ((2)1. The circles correspond to the two stellar 
models (e = and 0.8) at the point X in Fig.m 



nonrotating background quantities of NM and DM sep- 
arately. However, their definitions are meaningful only 
when the two fluids are non-interacting (as we assume 
in this work). In the general situation, the angular mo- 
mentum of each fluid would contain a contribution which 
depends on the coefficient A in Eq. ([1]) and the relative 
velocity J^nm — ^^dm [13 • The coefficient A vanishes only 
when the master function A is independent of the scalar 
product = —Uap". 

With the moments of inertia of NM and DM defined 
individually as above, we can further define the total mo- 
ment of inertia of DANS by / = /nm + /dm ■ In Fig. 21 we 
plot the scaled moment of inertia / against compactness 
M/R for the DANS sequences shown in Fig. [T] As in 
Fig. [l] the solid lines are sequences for different amount 
of DM specified by e. The case e = corresponds to 
ordinary neutron stars. The dashed line corresponds 
to Eq. ([2]). The three vertical lines (with arrows) at 
M/R — 0.05, 0.1 and 0.15 represent the range of val- 
ues of / obtained by the large set of EOS models which 
were used to obtain Eq. ^ . They can be regarded as the 
error bars of Eq. ([2]) at those values of M/R. The cir- 
cles in the figure correspond to the ordinary neutron star 
(e = 0) and DM dominated star (e = 0.8) at the point X 
in Fig. [TJ While the scaled moment of inertia of ordinary 
neutron stars can be modeled approximately by Eq. ([2]) , 
Fig. S] shows that / of DANS depends sensitively on the 
amount of DM. In particular, for the DM dominated se- 
quence e — 0.8, the value of / is significantly smaller 
than that allowed for ordinary neutron stars with the 
same compactness. This might lead to observational sig- 
natures of DM dominated compact stars. Further work 
is required to investigate the observational implications 
of DANS in details. 

So far we have focused our study on DM particle mass 
mx — 1 GeV. In view of the recent interest in DM candi- 
dates in the mass range of a few GeV, it is interesting to 
see how different DM particle masses in this range affect 
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FIG. 5: Maximum stable mass A/max is plotted against the 
DM particle mass mx for a fixed amount of DM specified by 
e = 0.1. 



Discussion and Conclusions. — In summary, we have 
studied the effects of a degenerate DM core formed by 
non-self-annihilating DM particles of mass ^ 1 GeV upon 
the structure of neutron stars. The structure of these 
DANS depends strongly on the size of the DM core. In 
particular, we found a new class of compact stars which 
are DM dominated — a NM core embedded in a ten- 
kilometer sized DM halo. The stability of these stars 
has been checked by performing a radial perturbation 
analysis. These DM dominated stars have rather differ- 
ent mass-radius relations and (scaled) moments of inertia 
comparing to ordinary neutron stars. A distinctive prop- 
erty of these stars is their small NM core radius of about a 
few km, from which thermal radiation could be observed. 
The detection of a compact star with a thermally radi- 
ating surface of such a small size could provide a strong 
evidence for their existence. 



our results. In Fig. [51 we plot the maximum stable mass 
-Mmax along a sequence of compact stars with e = 0.1 as 
a function of mx- Note that the sequence for mx = 1 
GeV is shown in Fig. [TJ For a given proportion of DM 
inside the stars, Fig. [5] shows that a higher DM particle 
mass in general leads to a smaller maximum stable mass. 
This can easily be understood by noting that, since the 
NM and DM are assumed to be non-interacting (except 
through gravity), the DM core is supported only by its 
own degenerate pressure. It is well known that the maxi- 
mum mass limit for a self-gravitating Fermi gas decreases 
as the particle mass increases. Hence, the onset of the 
collapse of a degenerate DM core is responsible for the 
dependence of M^ax on mx as seen in Fig. [5] It should 
also be noted that, while the pressure of NM within the 
DM core does not contribute to supporting the weight of 
the DM core, the mass of the NM fluid does enhance the 
collapse of the DM core. 



Could DANS be formed in the first place? To answer 
this question, one needs to consider the effects of DM 
on the stellar formation process. In fact, the heating ef- 
fects due to DM annihilation on stellar formation have 
been studied in recent years [22| . The result is the pre- 
diction of a new phase of stellar evolution during which 
a protostar is supported by DM heating. What if one 
replaces the annihilating DM model in 22\ by non-self- 
annihilating DM? How would such a DM core affect the 
stellar evolution? Would the DM core survive the su- 
pernova explosion of massive stars and form DANS as 
studied in this paper? These are challenging questions 
that deserve further investigation. 
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